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We develop an analytical method for calculating the geometric wakefield and impedances of an 
ultrarelativistic beam propagating on- and off-axis through an axially symmetric geometry with 
slowly varying circular cross-section, such as a transition. Unlike previous analytical methods, our 
■ approach permits detailed perturbative investigation of geometric wakefields, and detailed pertur- 

^SJ , bative investigation of impedance as a function of frequency. We compare the accuracy of the 

. results of our approach with numerical simulations performed using the code ECHO and determine 

' parameters in which there is good agreement with our asymptotic analysis. 

>• ! PACS numbers: 41.20.Jb, 41.60.-m 
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I. INTRODUCTION 



Vital considerations in modern accelerator and light source design include the influence of non-uniform metallic 
structures, such as a vacuum system or collimator jaws, on nearby charged particle beams. Rapid changes in the 
spatial profile of the structure tend to have undesirable consequences for a particle beam, such as inducing instabilities 
I '' and emittance growth; hence, designers often employ structures with cross-sections that gradually vary with distance. 
^ , For example, a series of gradually tapering metallic structures (collimators) may be used to strip unwanted particles 

■ from beams prior to the collision event (see, for example, [ij). 

^ , Optimal design of particle accelerator subsystems is commonly sought by direct numerical solution of Maxwell's 
^ ' equations, but accurate numerical computation of the electromagnetic fields near a short bunch in subsystems such as 
a post-linac collimator requires large computing resources In particular, numerically resolving a collimator with a 
gradually tapering geometry requires a mesh whose cells are much shorter than the length of the collimator. Although 
f~| ' one may employ windowing techniques to avoid calculating the fields throughout the entire structure at every time 
O-i step, such calculations frequently require intensive parallel computation [j, Similar problems are encountered 
when considering beam pipe transitions inside small-gap undulators in future fight sources. In such cases, analytical 
i expressions for the beam's behaviour are highly desirable. 
^ ' In practice, it is assumed that the beam does not deviate much from rectilinear motion parallel to the axis of the 
00 . metallic structure The beam's trajectory is often obtained as a perturbation due to the coupling impedance of 

■ the unperturbed beam and structure (early discussions of impedances in tapered structures were given by Yokoya Q 
] and Stupakov Q)- Although each frequency component of the unperturbed beam leads to its own impedance, it has 
, been noted that low-frequency transverse impedance is dominated by the zero-frequency component of the source 

■ and longitudinal impedance is directly proportional to frequency over a broad frequency range 7]. Recent work [I,!] 
, has been tailored to the above observations; it yields a frequency independent result for the transverse impedance 

and a longitudinal impedance directly proportional to frequency. Such methods do not permit detailed exploration 
of impedances as a function of frequency; the purpose of the following is to address this limitation. 

To overcome the above, a new scheme for obtaining analytical expressions of impedances in gradually tapering, 
axially symmetric structures was suggested by us in (lo| . In common with Q , our approach employs an expansion in 
a small parameter e characterising the gradually changing cross-sectional radius of the structure. However, we employ 
a decomposition of Maxwell's equations using auxiliary potentials that yields impedance as a series in frequency, 
without a priori assuming that the frequency is small. 

The following is an extensive investigation of the approach introduced in (lo| . We study the longitudinal and 
transverse wake potentials of a bunch travelling parallel to the axis of a perfectly conducting axially symmetric 
structure and compare the second, fourth and sixth order results in e with numerical data from ECHO [11]. We then 
present analytical expressions for longitudinal impedance; our method applied to a harmonic source reproduces the 
results given by Yokoya @ and Stupakov 0, 0] when working to second order in e. We also present corrections to 
the Yokoya-Stupakov results that arise due to higher order terms in e. Expressions for the longitudinal impedance 
up to fourth order in e have been given previously in |l5 | and we include them here for completeness. We then turn 
to a detailed exposition of the passage from Maxwell's equations to our iterative procedure for calculating auxiliary 
potentials order-by-order in e. A discussion of the difficulties encountered when auxiliary potentials are not used is 
given in Appendix IB] 



FIG. 1: A waveguide with a slowly varying profile 



II. OVERVIEW AND RESULTS 



A full account of our solution method is presented in later sections. This section focusses on a comparison of the 
results obtained using our asymptotic method and results obtained using the code ECHO. 

Our investigation concentrates on beams propagating through axially symmetric structures whose circular cross- 
sections gradually change along the axis of symmetry. Such structures are topologically equivalent to an infinite right 
circular cylinder and we call them waveguides (see Figure [T|). This article focusses entirely on the geometric wake 
of the beam in the inductive regime Here, the walls of the waveguide are assumed to be perfectly conducting; 
resistive effects will be considered elsewhere. 

We develop the electromagnetic field of the beam as an asymptotic expansion in a small parameter e characterising 
the gradually changing radius of a waveguide with a smoothly varying profile. In particular, the radial profile of the 
waveguide is specified as 

r = R{z) = R{ez) (1) 

where e ^ 1 and (r, 0, z) is a cylindrical coordinate system with z a Cartesian coordinate parallel to the waveguide's 
axis of symmetry (r = 0). 

The two radial profiles explored here numerically are 

r = R{z) = 20-18 exp[-zV(8 x lO^P)] (2) 

with e = l/\/8 X 10^ and 

r = R{z) = 20-18 sech (O.Olz/0 (3) 

where e = 0.01. The coordinates z, r and the constant I have dimensions of length; in the following I = 1mm, and z, r 
are measured in mm. The profiles ([21 El have a ratio of gap to beam pipe radius and average taper gradient similar 
to that proposed for the beam delivery system of the next generation of lepton colliders [l3] . 

Wake potentials are calculated by integrating components of the electromagnetic field along a straight line parallel 
to the axis of the waveguide. We consider fields excited by a narrow Gaussian bunch of effective length a propagating 
at the speed of light parallel to the waveguide's axis, offset by a distance tq from the waveguide's axis. The charge 
density p of the bunch is 

q^e 2cr^ 

Pro,Oo i.'T, ^, u) = " COS 6*0) (5 (y - ro sin 6*0) (4) 

(TV 27r 

where u = z — ct^ x = rcos9, y = rsinO and go is the total charge of the bunch. The longitudinal wake potential M^" 
is given by 

1 f°° 

WHr,9,u) = --l i?4,_dz (5) 



00 



where 

E,l^^^_^=E,_[r,e,z,t=^j (6) 

and Ez is the longitudinal (i.e. z) component of the electric field. The transverse wake potential Wj^ may be obtained 
from MK'I using the Panofsky-Wenzel relation p^ . 

Sample results for ly" are shown in Figures [2l [SI The "leading order" curves result from electromagnetic fields 
calculated up to 2nd order in e, while the "4th order" and "6th order" curves arise from including 4th order and 6th 
order corrections in e respectively. It may be shown that the 3rd, 5th and 7th order contributions to VF" vanish. 
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Wake potentials can be calculated from impedances using Fourier methods [l^ and the leading order results may be 
recovered using well-known expressions for impedance originally due to Yokoya Corrections to the leading order 
transverse impedance that properly accommodate the waveguide's taper but neglect the frequency dependence of the 
source may be found in Q; we expect the frequency dependences of impedances to be more accurately represented 
using our method. 

Data from ECHO was obtained for a sequence of mesh cell densities (10, 15, 20 cells over bunch length cr; we judged 
the numerical errors to be small based on the small differences between the results for 10, 15 and 20 cells per cr). The 
monopole and dipole contributions to ^ are the zeroth and first order terms in an expansion of ^ with respect to 
To; for direct comparison with the ECHO data, only the wake potentials due the monopole and dipole contributions 
are calculated. 

Figures (|3117]) show the magnitude and position of the maxima of longitudinal wake potentials, the magnitude of 
the maxima of the radial component of transverse wake potentials and the momentum kicks induced by transverse 
wakes in the exponential geometry. Figures (|8l fTT|) show the corresponding quantities for the sech profile. The higher 
order contributions reduce in significance as the effective bunch length a increases and we observe that the 4th and 
6th order results agree well with the leading order predictions for bunch lengths longer than ^ 2mm. For a less than 
~ 1.5mm the data shows that, in a number of cases, the maximum values of the wake potentials and momentum kicks 
obtained using ECHO generally agree better with the 4th order results than with the leading order (Yokoya-Stupakov) 
predictions over the parameter ranges considered here. The sample plots of M^" shown in Figures [H [3] demonstrate 
this for a bunch length of 1.2mm; for such bunch lengths the reliability of the 6th order correction is questionable. 
Although the agreement with the numerical results is impressive, we warn the reader that the radius of convergence 
of our perturbative expansion is unknown. 



A. Impedance of an on-axis beam 



In this section we present impedance formulae due to the on-axis harmonic charge density 

g^{r,e,u) ^ —e'^S{x)S{y) (7) 
c 

where /^^ defines the harmonic current component of the beam. The longitudinal impedance Z^^^lu) for angular 
frequency uj is 

ZH^) = -Y (8) 

J —oo 

where is the longitudinal component of the electric field generated by the source ^ . The approach detailed in 
subsequent sections leads to the following : 

^on-axis ('^) "^Ion-axis ^ on-axis ^ "^4 on-axis ^ on-axis • ■ ■ ('^) 

where z\^^_^^:^^ + ^J^n.^xis ^^'^ Yokoya-Stupakov longitudinal geometric impedance 



^2onaxis = -:^ ^ ^ (11) 

ion-axis 47r£(,c2 J ^ ' 

with B! — dR/dz, Ri — R{oo) and i?2 — R{—oo). Equations pUl [TT|) follow, respectively, from 1st and 2nd order 
terms in e in an asymptotic series for E^. The corrections 

4^socZi„.,x,, = £ 1^ |5i?'^ + 3 iRR"f - 2^ (i?^ R'f |d. (12) 



and 



— (i?3 R' R"f + J^i?7 ^"3 „ Uz (13) 

160^ ^ 1920 1920^ ' ' ^ ' 
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FIG. 2; Sample longitudinal wake potential W" versus a = z — ct for an on-axis source (length 1.2mm) propagating through 
the exponential geometry ([2}. 

arise from 4th and 6th order terms in the asymptotic series for Ez- The imaginary part of is invariant under 
reversal of the direction of the beam, i.e. Im(Zll ) is invariant under the replacement R{z) i-^ R(—z). Thus, each term 
in the integrand in Im(Zll) only contains an even number of derivatives of R and arises from an even power of e in 
the asymptotic series for Ez- It follows that the odd order contributions to Im(Zll) vanish. 

B. Impedance of an ofF-axis beam 

This section focusses on longitudinal impedance due to the off-axis harmonic line charge density 



guj,ro.ea {r, 9, u) = — e"^"(5 (x - tq cos 6*0) S {y - ro sin 6*0) 



(14) 



where rg 7^ 0. The corresponding transverse impedance Z_l may be obtained using the Panofsky-Wenzel relation [l3 |. 
Fourier expansion of the source in 9 leads to 



Zl +Z!j+Zl+ Zl 



(15) 



where ^on_axis (^^^ ®) monopole contribution to and the contributions arising from 1st, 2nd, 4th and 

6th-order terms in E, are 
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FIG. 3: Sample longitudinal wake potential W" versus u = z — ct for an on-axis source (length 1.2mm) propagating through 
the hyperbolic secant geometry (jS)). 
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FIG. 4: The maximum in u oiW" versus effective bunch length a for a monopole source propagating through the exponential 
geometry ((2|. 
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4' = 7^ i: 3 ( i " i ) r-'r^ cos m(0 - ^o), 



47r£oC ^ m V i?? 

m— 1 \ J- 



m— 1 

oo 



II ^-^ 2ia; r^'r™ cos TO (0 — 6*0) 
AtteqcZI = 2^ 



i?"^ 2TO^ + 6m+l i? 



//2 



c to(to+1)(to + 2) 1 y_oo V-R^™ 3 
uj^ 3TO2 + 8m + 6 (Am -2, R"^ R"^ 



:2 to(to+ 1)2(to + 3) V 3 i?2m-2 ^2m-4 



dz 



(16) 
(17) 



(18) 



I X ^ iuj r^r™ cos to(6' — 6*0 ) 
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respectively, where to is the multipole index. Since lm{Z") is invariant under the replacement R{z) 1— > R{—z), i.e. 
invariant under reversal of the direction of the beam, the odd order terms in the asymptotic expansion of Im(Zll) in 
e vanish. 

The longitudinal impedance of a dipole source is obtained by setting ni = 1 in ([TBI [TTI [T8l [T9|) : 



47reocZ|' = 2 ( i- - i- ) rgr cos(0 - ^o), 



47reocZ2 — — 2i— ror cos(6' — 0o) / ^^dz. 
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FIG. 5: The position of the maximum in u of W" versus effective bunch length a for a monopole source propagating through 
the exponential geometry ([2}. 



Use of the Panofsky-Wenzel relation [Tj| and ((201 [21] [22l [23| leads to the following contributions to the transverse 
dipole impedance 
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Z^(cj) = + Z.^ + zi + zi 



(27) 



(28) 



The imaginary part of the transverse impedance at zero frequency lim^j^o Ini[Z''"(w)] was previously obtained in 
As far as we are aware, the above expressions for transverse impedance as a function of frequency are new. 



III. DERIVATION OF RESULTS 



The remainder of the paper focusses on the method used to derive the expressions presented above. 
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FIG. 6: The maximum in u of the radial component of Wx versus effective bunch length a for a dipole source propagating 
through the exponential geometry ([2]). 
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FIG. 7: The transverse momentum kick versus effective bunch length a for a dipole source propagating through the exponential 
geometry 



9 



Maximum Wake vs Bunch Length 




0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 

Sigma (m) 



Leading Order — • — 4th Order — • — 6th Order 
ECHO 10 cells/Sigma — • — ECHO 20 cells/sigma 



FIG. 8: The maximum in u of W" versus eflective bunch length a for a monopole source propagating through the hyperboUc 
secant geometry Q. 
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FIG. 9: The position of the maximum in u of M^" versus effective bunch length o for a monopole source propagating through 
the hyperbolic secant geometry (|3]). 
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FIG. 10: The maximum in u of the radial component of Wx versus effective bunch length a for a dipole source propagating 
through the hyperbolic secant geometry 
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FIG. 11: The transverse momentum kick versus effective bunch length a for a dipole source propagating through the hyperbolic 
secant geometry (|3]). 
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A. Decomposition of Maxwell's equations using auxiliary potentials 

1. 2+2 split of Maxwell's equations 

The approach we use is based on a 2+2 split of Maxwell's equations adapted to the waveguide. We employ cylindrical 
polar co-ordinates {t, r, 6, z) on Minkowski spacetime where r = is the waveguide's axis, 9 is the azimuthal angle 
around the waveguide's axis and z is a co-ordinate parallel to the waveguide's axis. We employ exterior differential 
calculus to analyse Maxwell's equations (see, for example, [H, [l^) as it affords a concise language for exploiting 
different coordinate systems and for developing our new auxiliary potential method discussed in the following. 

Let {A4,g) be Minkowski spacetime where g is the metric 

g = -c^dt «) dt + dz «) dz + dr (K) dr + r'^dO (g) d9 (29) 
in {t, r, 6, z) co-ordinates. We choose the volume 4-form ^1 as 

*1 = cdtAdzA#il (30) 

where the 2-form #_l1, 

#_l1 = r^dr Ad6', (31) 

and the metric gx, 

g^=dr®dr + r^d9 (g) d0, (32) 

are adapted to constant {t, z) surfaces inside the waveguide. The waveguide surface in is the 3-dimensional 
hypersurface {r = R{z) \ — oo < z < oo} where R{z) is the radius of the cross-section at z with unit normal dz. 

The source moves in the positive z-direction near the speed of light in the laboratory frame, and its 4- velocity field 
V is well-approximated by the null vector : 

V = -dt + dz. (33) 

c 

The metric isomorphism between the spaces of vectors and co- vectors on spacetime is denoted as follows: the co- vector 
field X associated with the vector field X is defined such that X{Y) — g{X,Y) for all vector fields Y and the vector 
field 5 associated with the co- vector field a is defined such that /3(5) — g^'^{(3,a) for all co- vector fields /3. Hence, 
the 1-form V associated with V is 

V = dz- cdt. (34) 
It is expedient to transform to a new co-ordinate chart: 

u = z — ct, C, — z (35) 

with 

g = dC, ® du + du ® dC, — du ® du + g (36) 
★ 1 = dCAduA#il. (37) 

The coordinates psp play a useful role because the source is naturally expressed in terms of u and the boundary of 
the waveguide is described using Q. It follows that 

V = du, (38) 

*i/ = duA#j.l. (39) 

Here, the Hodge map 7k- on satisfies 

★ (j/ A X) = ix * (40) 

where v is an arbitrary form and X is an arbitrary vector field and ix is the interior contraction with X. The 
transverse Hodge map is defined on the transverse, cross-sectional subspace to satisfy 

#x(/i±Ay)=iy#xM± (41) 
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where is a form on the domain spanned by {dr, d9} and F is a vector field in the span of {dr, dg}. 
The Maxwell equations on space-time are 

di^ = (42) 
eodirF = -gi^V (43) 
where gV is the electric 4-current of the beam and d is the exterior derivative on forms. Charge conservation 

d{g*V)^0 (44) 

yields 

d^g^O. (45) 

Without loss of generality, the electromagnetic 2-form F is decomposed as 

F = $dCAdii + dwAa + dCA/3 + *#±l (46) 

where a and /3 lie in the subspace of 1-forms spanned by {dr, dd} (transverse subspace). The component $ of is 
the longitudinal component of the electric field in (r, 0, u, C) coordinates: 

E. (^r, 0, z = C, t = ^(C - ")) = 9, u, C). (47) 

(See, for example, [l3] for a detailed discussion of the relationship between F and electric and magnetic fields in vector 
notation). It follows that the Maxwell equations (|42)) and (|43|) can be decomposed as 

d^<i> + d^a - duP ^ 0, (48) 

ac*#il-d^/3 = 0, (49) 

a„*#il-d^a = 0, (50) 

ac#i(a + /3)+5„#i/3-cU* = 0, (51) 

di#i(a + /?) - a„$#il = -pir,9,u)#^l, (52) 

9c$#±l + di#i/3 = (53) 

where d^ is the exterior derivative acting in the transverse subspace, and p{r,9,u) = j^g{r,9,u), the functional 
dependence of the charge density being written explicitly to emphasize its independence of ^ due to (j45p . 

2. Auxiliary potentials 

In this section we introduce a set of six 0- forms (auxiliary potentials). This anticipates the imposition of the 
perturbation scheme in Section llll Dl they will be shown to satisfy tractable second order PDEs in gradually tapering 
waveguides, whereas working directly with the conventional electromagnetic potentials leads to third order PDEs. 
Appendix [B] contains a brief discussion of this issue. 

After rewriting the Maxwell equations in terms of auxiliary potentials, and adopting perfectly conducting boundary 
conditions, we will show that asymptotic expansions of the auxiliary potentials yield 2-dimensional Poisson and 
Laplace equations at each order in the expansion parameter e. 

In Appendix \^ we demonstrate that the transverse 1-forms a and /3 can be replaced by the 0-form pairs {A, a} 
and {B, b}, viz. 

a = d^^ + #_Ld_La, (54) 
f3 = d^B + #^d^b (55) 

where A and B obey Dirichlet boundary conditions on the waveguide boundary, i.e. A = B = Oa,tr = RiC)- 
It follows that (since d^ = 0) the Maxwell equations (|48ll53p partially decouple to give 

di_ ($ + dcA - duB) + #^di_ idea - dub) = 0, (56) 

+ 6j_d^b = 0, (57) 

9„* + (5j_dia = 0, (58) 

d± {dciA + B)+ duB) + #_Ld_L (* + d^ia + b) + d^b) ^ 0, (59) 

6±d±{A + B) + du^^ p{r,e,u), (60) 

9f $ - 6±d±B = 0, (61) 
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where 

S^fi^ = (-l)P#lMi#iM± (62) 

is the transverse co- derivative of a transverse (i.e. in the span of {dr,d9}) p-form fij^ and #J^#_l = #±#1^ = 1- 
Motivated by the appearance of equations ([55)1 and introduce 0-forms {7i*, H''^, , H'', W, X} that satisfy 

duH'^ = $ + - duB, (63) 

duH"^ = dca - (64) 

duH^ ^d(;{A + B)+duB, (65) 

7^'' = ^- + 9^(0 + 6) + 9„6, (66) 

duW^A + B, (67) 

5„X = a. (68) 

Equations (|63ll68p do not specify {H"^ ,W , X} uniquely, as they all appear under a w-derivative. Additional 

"gauge fixing" conditions will be supplied in the following. 
Substituting ((67|) into ([65]) and integrating with respect to u gives 

B = H^ ~d^W + X{r,e,C) (69) 
where X{r,9,() is an arbitrary 0-form of integration and, hence, from [S71 [S^ . 

A = d^w + dcW -n^ - X{r,eX), (70) 

$ = + + d^n^ + 9c^(r, 0, C) - 2dl^W - d^^W. (71) 
Substituting into and integrating with respect to u gives 

b^dcX-H^-j{r,e,C) (72) 
where 7(r, 9, C) is a 0-form of integration. Substituting (|68p and ([72)1 into ([66)1 gives 

* = -2d^^X - dl^X + dcU"^ + SuW'^ + 9c7(^, 6, () + H^ (73) 
We can simplify the above by introducing 0-forms A(r, 0, C,) and r(r, 0, 

9cA(r,0,C) = A(r,0,C), (74) 

9cr(r,0,C) -7('^,^,C) (75) 

and 0-forms W , X, 

W ~ A{r,0,C), (76) 

X = X-Tir.,9,C). (77) 

Substituting the above into (|68ll73p gives 

B = n^ - dcW, (78) 

A^duW + dQW -n^, (79) 

$ = + a^W^ + c»cW^ - 2dl^W - dliW, (80) 

a = (81) 

6 = - W^, (82) 

* = -2a|„A: - dl^x + dt^n.'f + duH'^ + n\ (83) 

Having derived a mapping from the auxiliary potentials {H^ ,W,Ti'^ , X,Ti'^} to {B, ^, a, 6, we now rewrite 
equations (|56ll6ip in terms of {H^, VF, H*, X, H*^}. Equations ^ and JH]) become 

a„ (d^H* + #±diW'^) = 0, (84) 

dudiJi^ + #±d±n'' = 0. (85) 
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Equation (IM)) is equivalent to 

diH-^ = #idiH* + K(r, 0, C) (86) 
where K(r,9,Q is an integration 1-form. Decomposing k as (see Appendix 

^(r, C) = #±di «i (r, 0, C) - d^K2 (r, 0, () 

and using ([5^ yields 

d± (Ti'^ + K2(r, 0) = #±di (W* + K, (r, 0, 0) . (87) 

It can be seen from (I78II83|) that the^O-forms A, a, B, b, "if (and hence our electromagnetic 2-form F) are invariant 
under the "gauge" transformations [22] 

n"" ^n"" + dcQ{r,ex), (88) 

W ^W + Q{r,e,C)+w{r,9), (89) 
W*^W*+n(r,6i,C), (90) 

rc^ ^rc^ + dcM{r,e,c), (9i) 

X + M{r, 0, C) + x{r, 9) (92) 

where Q, w, n, M and x are arbitrary functions of the indicated variables. In what follows, these "gauge functions" 
will be chosen to reduce the Maxwell equations to a form amenable to the gradual taper approximation. 

Without imposing any restrictions on the electromagnetic 2-form F, we choose {A/, n} so that K2} satisfy the 
"gauge" conditions 

5cM(r,f?,C) = -A^2(r,^^,C) (93) 
n(r,0,C) = -Ki(r,^?,C). (94) 

Equation ([87)1 becomes 

= #±diW*. (95) 
Applying the transverse exterior derivative d^ and co-derivative (5^ to (|95p gives 

JidxTi* = ^id^W'^ = 0. (96) 

Similarly, equation ([55]) implies 

JidiW^ =a(r,^,C) (97) 
where (T(r, 9, Q is a 0-form of integration and so, by choosing Q{r, 9, ^) such that d(^6±d±Q{r, 9, () — a{r, 9, Q, we find 

dj_di_n^ = 0. (98) 

Thus, {7i*, W', , H^} may be chosen to be harmonic with respect to the 2-dimensional transverse Laplacian. 
Now turning to ([571 [55)) . we note that the Maxwell equations (|58l) and ([5T|) give 

du {6±d^X + *) = (99) 

so that 

S^d^X + ^ = ,y{r,9,C) (100) 
where ;/(r, 9, C) is a 0-form of integration. As T-C^ is harmonic (see the Maxwell equations l|57p and ([5^ give 

5c (^±dj.A: + f) = (101) 

and it follows 

^id^AT + l' = /^(r,6',w) (102) 
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where /i(r, 6, u) is an integration 0-form. The left hand sides of (|100p and (|102l) are equal, so 



(103) 



Hence, the integration 0-forms ^(r^9,u) and i>{r,9,C,) are constant with respect to u, C- Eliminating 5* from pOOp 
using yields 



(104) 



where fi{r,6) = fi{r,6,u) = i'{r,d,Q- The function x{r,d) has not yet been used in the transformation ([92]) . We can 
eliminate /t(r, 0) without affecting the fields by setting 



Thus, we conclude that (jST)) and ((58)) may be written as the single equation 



S^d^X - 2dl^X - dlc^X + dc^H'^ + duH^ + - 0. 



(105) 



(106) 



We now consider the final two Maxwell equations (pDl . Using ([751 [7^ , the Maxwell equation ([50)1 may be written 
as 



du[5^d^_W + ^)^ p[r,e,u). 

Introduce P{r, 9, u) where 

duP{r,e,u) = p{r,e,u). 

Equation (I107p can thus be integrated to give 

6±d±W + $ = P(r, 6, u) + cr(r, 6, () 
where <;(r, 9, () is a 0-form of integration. Now, using ([78]) . and noting that is harmonic, ([61]) gives 

9c((5±dj.VK + $) =0 

which implies 

S±d±W + $ = x(r,d,u) 

where x('', ^','") is a 0-form of integration. Equating the right hand sides of (|109p and piip gives 

Thus, <;(r, 9, C) is constant with respect to ( and eliminating $ from (|109p using ()80p yields 

(5idxW^ - 2^2^ _ ^2^1^ + duH^ + d^n^ = P(r, 9, u) + <f(r, 0) 

where <f(r, 61) = <j(r, 9, (). Using the transformation ([89]) we may set 

5^d^w{r,9) = <;{r,9) 

and conclude that, without loss of generality in F, one may choose a "gauge" in which 



S^di_W - 2dLW - d?.W + duH"^ + duH^ + d<^n^ = P{r, 9, u) 



(107) 
(108) 
(109) 
(110) 
(111) 
(112) 
(113) 



(114) 



is satisfied. 
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3. Summary 

We seek appropriately bounded solutions to the following system of equations for our auxiliary potentials W, X, 
n^, n\ and W^: 

5i_dil-L^ = 0, (115) 
Ai^v!" - #±di , (116) 

d^n^ = #idiH*, (117) 

5^di_W - 2dlcW - dl^W + + duU^ + dcU^ = P(r, 9, u), (118) 
S^d^X - 2dl^X - dl^X + di^H.'f + duV.'^ +0!' = Q. (119) 

Notice that (fTTC)) and pTZ)) imply that S^d^Ti}' = S^d^K^ = S^d^H"^ = 0. The source term P(r, 6*, u) in equation 
(fTT8| satisfies equation (fTOSll : 

duP{r,e,u) = p{r,e,u) 

and the electromagnetic 2-forni F is 

F = $dC A du + du A (d_L^ + #±d_La) + dC A (d_LB + #±d_L6) + *#_l1 (120) 

where 

^ = 9„M/ + acM^-H^, (121) 
B = n^ ~ dQ,W (122) 

$ - + duH^ + 9c^^ - 29^c^ ~ ^CC^' (123) 

a = duX, (124) 

b^d^X-n'^, (125) 

* = - (2a^^X + 9^^X) + d^H"^ + duH''' + (126) 
Equations (jllSp and pi9p also give alternative expressions for $ and 

$ = P(r, e*, ?i) - (J^d^M^, (127) 
* = -(5_Ld^X (128) 

and F may be written 

F = + + dcH^ - 2dl(W - dl(W) dc a du 

+ dw A [di (a„T4^ + 9^T4^ - n^) + a„#id_LX] 

+ dc A [di (7i^ - a^vK) + #idi {dc^x - 

+ (9cH'^ + dnU"^ +1^^- 2dl^X ~ dl^X) #^1 (129) 

At first glance, it may seem that (|115mi9|) is no more amenable to analysis than the original Maxwell equations. 
However, we show in the following that these equations reduce to 2nd order PDEs that are straightforward to solve 
when the auxiliary potentials are slowly varying with Q. 



B. Boundary conditions on the auxiliary potentials 

The boundary of the waveguide is the surface 

r-i?(C)==0. (130) 
Given a function (/)(r, 0, (, u) such that (/)(i?(C), 9, (, u) = 0, it follows immediately that 

^0(i?(C),0,C,") =0, (131) 
^0(i?(C),e,C,") =0, (132) 

^</)(i?(C), 0, C, ") = + ^'(^)|:) '^(^(^)' ^' ^' = ^- (133) 
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where R'{C) = dR/d(. 

Orthogonal decomposition (see Appendix |^ of the 1-forms a and /3 employs potentials A and B that satisfy 
Dirichlet boundary conditions on r = R{C)- Hence, on the boundary (|12ip and (|122p give 



- d(;W = 



=i?(C): 



duW = 0\r^ii(^Q. 



(134) 
(135) 



Equation (|134|) serves as a boundary condition on , while ()135p is satisfied by imposing Dirichlet boundary 
conditions on W: 



(136) 



The usual conditions on the electric and magnetic fields at a perfectly conducting surface 5 = in spacetime may 
be written 



d5Ai^ = 0|5=o- 

Using S = r — R{() we can decompose (|137l) as 

^dr-R'{Ch = 0\r=R(c), 
drA/3 + *i?'(C)#±l-0|,^fl(C). 

With a given by ([5i)) . the d6 component of (|138p is 

deA + rdra = 0\r=R{Q. 

Since A satisfies Dirichlet boundary conditions at r = R{C) it follows 

rdra = ^ d^^X = 0\r=R(Q 
and this can be satisfied by imposing the boundary condition [23l | on X: 



The dr component of (|138p gives 



$ + i?'(C) -9ea-c),A = 0|, 



--RiO 



and using ([^ I124p to eliminate ^ and a, (|143p can be rewritten 

1 



where 



duB = 0|,^«(^), d^A + i?'(C)a,^ = Q\r=R[o 



have been used, which follow from p321 1133]) and the Dirichlet conditions A = B = Q\r=R{()- 
Equation pl4|) can be satisfied by imposing the following boundary condition on Ti. . 



(137) 



(138) 
(139) 



(140) 



(141) 



(142) 



(143) 



(144) 



(145) 



(146) 



We will now argue that (|139p is automatically satisfied without further conditions on {W, X^H ,H }. Since B = 
0\r=R(Q, it follows dgB = 0|r=_R((;) and (|139p becomes 



drb + R'iC)-^ = 0\r=Ric)- 



(147) 
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Using (|125p and (|128p to eliminate b and and noting 



equation (|147p becomes 
However, equation pi7p implies 



(148) 



and hence, from (I146p . 



d^n"^ = — dgH" 
r 



Furthermore, using (|1421 1133P with 4> — rdrX gives 



dLX + R'iQ-dr {rdrX) = 0|^. 



-RiO 



(149) 



(150) 



(151) 



hence (|148p is trivially satisfied using (|150( I15ip . Thus, (|139p is satisfied without further conditions on 

{w,x,n^,n'^}. 

In summary, the perfectly conducting boundary conditions on the electromagnetic field &i r — R{C) are satisfied by 
imposing the following boundary conditions on {W, X, , 7Y*}: 



W^O\r=R(Q, drX = 0\r^R(^c:), 



-R'iO-90X\r=R(Q. 



(152) 



C. Untapered Waveguide 

Our method for determining {W,X ,Ti.^ ,H'^} is based on asymptotic series for fields modelled on those in the 
untapered waveguide. 

When the waveguide is a regular cylinder with zero taper (i?'(C) — 0), the Dirichlet boundary condition on W 
yields d^W = 0\r=R (see (jlSSp ). Thus, ()152p immediately gives 

^n"^ ^0\r^R. (153) 

The auxiliary potentials {7i^,7i*} satisfy the transverse Laplace equation and {7i^,7i*} vanish at r = R{C)- Thus, 
it follows that 

/ ds_n^ A #^d^H^ = [ H^i^^d^n^ 
Jv Jav 

- [ n^d±#±d±n^ 

Jv 

= (154) 

where "D is the cross-sectional disc domain, and is its boundary. This implies that is constant on V, and 
therefore vanishes. A similar argument applies to 7i*. Furthermore, using ()1161 fTT7|) it follows that {T-C^,H'^} are 
constant on V and also may be chosen to vanish. Hence, (|118p and (I119P become 

S^d^W -2d^^W - d^^W ^ P{r,e,u) (155) 
6±djLX - 2d^^X - d^^X = (156) 



19 



and are to be solved subject to the Dirichlet boundary condition at r = RforW and the Neumann boundary condition 
at r = i? for X. For the present purposes, we are not interested in the source- free modes of the waveguide, so we set 
X = 0. Furthermore, charge conservation impUes that p is independent of ( and we expect the fields to be independent 
ofC 

For a source given by 

p = f{u)pir,0) (157) 

the 0-form P may be written 

P(r,9,u) — ^(u)p where -^S(w) — f{u) 

au 

and a solution to (|155p is 

W = E{u)W{r, 9) where S±d±W{r, 9) = p{r, 9) 

for W obeying the Dirichlet boundary condition at r = i?. The electromagnetic 2-form F is 

F = f{u)du^AlW (158) 

which is compatible with the usual assumption that the electroma gnet ic field vanishes ahead of and behind an ultra- 
relativistic source with compact support in u (see, for example, (sl, Il8l|l. 



D. Waveguide with a gradually changing radius 



In this section we develop asymptotic expansions of solutions to the preceding equations in an axially symmetric 
waveguide whose cross-section is a slowly-varying function of z = C. 

A waveguide is considered to be "slowly varying" in C, if ^(C) = R{^C) where e > is a small dimensionless 
parameter. Hence, the waveguide boundary is described by 

r-i?(eC) = 0. (159) 

Introduce a "slow" longitudinal co-ordinate s, which is defined by 

s = eC- (160) 

Rewrite all potentials as functions of s, using the notation 

x{r,9,eC,u) = x{r,9,s,u) (161) 
dcx{r,9,eC,u) = edsx{r,9,s,u). (162) 

A prime on a function accented with a caron denotes differentiation with respect to s. For example 

^m) = = ^^'(*)- (163) 



Our approximation scheme follows by writing 



ji=0 

where 

X e {A, a, B, b, !>, W, X, H^,H'', H*, H^} . (165) 

Substituting such series into (|115mi9|) and the boundary conditions (|152p we obtain PDEs at each order n. 

Since £?^x — equations (|115mi9p yield a set of transverse Laplace and transverse Poisson equations at every 
order n, and the boundary conditions on "H^ and 7i* now depend on {n — l)-order potentials. This leads to a 
straightforward procedure for calculating the potentials order- by-order. For each n we 
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1 . Calculate the harmonic potential Ti.^ by solving the 2-dimensional Laplace equation 

S^di_n^ = (166) 

subject to the boundary condition [24l| 

n^^W:,_, at r = i?(s). (167) 

2. Solve m 

dxH^^ = du#±di_n^ (168) 

for nl 

3. Calculate the harmonic potential 7i* by solving the 2-dimensional Laplace equation 

d^d^^nt = (169) 

subject to the boundary condition 

Ht = -R'{s)^dgXr,^i at r = i?(s). (170) 

4. Solve [H 

diW^ = #idi7i* (171) 

for w;?. 

5. Calculate the potential Wn by solving the 2-dimensional Poisson equation 

d^di_Wn = Wl_^ + 2a„VK;_i - duUl - djil - Ul_^ + P„ (172) 
subject to the Dirichlet boundary condition 

= at r = i?(s) (173) 

where 

p„(.,M)=|f^'''") j'^'^":^ (174) 

10 for n 7^ 

and duP{r, 0, u) — p{r, 0, u). 

6. Calculate the potential X„ by solving the 2-dimensional Poisson equation 

S^d^Xr. = K-2 + -'ft- duUl - 'Hl-i (175) 

subject to the Neumann boundary condition 

drXn = Q dXr^R{s). (176) 

We construct A„, a„, 6„, $„ and from 

A„ = a„w'„-HW^;_i-7i:f, (177) 

= Tif - W^„-i', (178) 
= -W^"-2 - '^duW'^-i + a^H* -f a.Hf + Tif./ (179) 
= P„(r,0,u)-(5xdiVK„, (180) 







(181) 


K 




(182) 






(183) 




= -5xd^X„. 


(184) 



Finally, we perform the summations to the required order and change variable from s to 

x{r,0,s,u)^xir,0,eC,u), (185) 
x\r,0,s,u)^e-'dcxir,e,eC,u). (186) 

Equation (jl20p then yields an asymptotic approximation for the electromagnetic 2-form F with e <C 1. 
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E. Sources cixially symmetric with respect to the waveguide 

When the source has rotational symmetry about the axis of the waveguide, the fields sought are independent of 9 
and determining auxiliary potentials independent of 6 is straightforward. 

For each order n: 

1. is constant in the transverse waveguide cross-section, and is equal to the boundary value of 



'^n — ^n-l\r=R{s)- 

2. The other harmonic functions, and X„ are zero, 

n= rLn = rty = /t„ = U 

and hence 



= &n = *n = 0. 



3. Wn is calculated via the integral 



lR(s) 



f)„ridri 



dr2 
r-2 



where 



4. Finally, B„, $„} are 



An = d^Wn + Wn-l -n^, 

^n = -w':,_^-2d^w'^_, + duK + K-i' 

= Q -\- P 

— ^^n ^ n' 



For example. 



$0 =0, 
=0, 



W 
H 

w 



= 



277 

277 ^(s)' 

f{u) R'(s) 



Inr-ln^(s)) , 



1 = — 



2^ R{s) ' 



Stt i?(s) 



(187) 

(188) 
(189) 

(190) 



(191) 
(192) 

(193) 
(194) 



(195) 
(196) 

(197) 
(198) 

(199) 

(200) 



- 

$2 = 



47r 

f{u) ifR'is) 



W 



^ 4^ \\r{s) 
/'(«) 



647r 



{R{s 



Ris) 



{3R{& 



')+4R'{sf^ 



(201) 
(202) 

(203) 
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where f'{u) = df{u)/du. 

The longitudinal impedance is 



1 r°o 

-'ti' n — oo 



and employing the above iterative procedure with the harmonic profile 



f{u) = — e 



leads to 



$0 =0, 

$2 



27reoc 



{HR{s)]y, 



e*"" j/i?'(s) 



327reoC'^ 

Hence, introducing -^Ion-axis where 

z" 



2 ^2) 



327r£oc'^ 



n on-axis 



(204) 

(205) 

(206) 
(207) 

(208) 



(209) 



(210) 



it follows dMl EOS) yield m E]). Clearly, Zjl' 



i'on-axis " ^ siuCC WC choOSC R' (oo) = R'{ — (X)) = 0. 



on-axis 



= and the third order contribution (|209p to $ leads to 



Although expressions for $„ rapidly increase in complexity with increasing n, they follow directly using the above 
iterative procedure and are straightforward to generate using a computer algebra package such as Maple [T^. 
It may be shown that {Z, 

{M' on-axis: 4'on-axis> ^io„-axis> 4'o„-axis} ^^6 givCU in ScctionHIl 



Son axis' ^7 on axis} vauisli bccause -R'(oo) = R'{—oo) = 0. Exphcit expressions for 

■II 



F. Source Offset from the Waveguide's Central Axis 



The general scheme developed in Section IIII Dl simplifies when used to calculate the geometric impedance of an 
infinitesimally thin beam offset by a displacement from the central axis r = 0: 



p{r,9,u) = f{u)S {x - ro cos 6*0)) (5 (y - ro sin 6*0) ■ 



(211) 



We first calculate the zero order fields. As before, since and Tig are harmonic and vanish on the waveguide 
boundary, it follows = Ti* — and we choose TYq = = 0. Equation (|172p gives 



^^d^T^o = ^(u)^ (x - ro cos 6*0)) S {y - tq sin 6*0) 
and a solution to (|212p which vanishes at r = R{s) is 



^ + R{sf ~ 2rro cos(6' -60)] - In (r^ + - 2rro cos{e - 6*0)) 



(212) 



(213) 
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It follows 



47r 



2 In i?(s) + In 



1 - ^e-^ 



R{sy 



In 



1 - 



where . . . indicates terms independent of s. For vq < R{s) we have 



and using 



i?(s)2 



< 1 



(214) 



ln(l - a;) = - 



rr. 



for Ixl < 1 



it follows 



Wp 



27r 



In i?(s) 



E 



1 



ror 



'-^^ m \R(s)^ 

where . . . indicates terms independent of s. Differentiating Wq with respect to s yields 

ror 



_E{u)fR'{s) ^^R'js) ^ 



costo(0 — 6*0) 



(215) 



(216) 



The first term in (j216[) is the contribution from the monopole (m = 0) term in the source and has the same form as 
Wq due to an on-axis thin beam. Thus, the monopole contribution has already been discussed in Section III Al and 
here we concentrate on the terms that arise for m > 0: 



W' 



0,m 



2m+l ^ ^ ^' 



TT i?(s)2™+ 



(217) 



The final potentials and impedances can then be evaluated by summing over m. A commonly used approximation for 
fields with r, tq ^ R{s) is to truncate the multipole series at the dipole contribution (m = 1); see, for example, 0]. For 
each order n > 0, the 6-step procedure of Section llll Dh nust be followed to obtain the auxiliary potentials. In practice, 
this procedure is simple to implement for an infinitesimally thin beam and the main features of the calculation are 
summarized below. 

Steps 1 and 3 are satisfied by and Ti* of the form /i^^*(s, u)r™ cos to(0 — ^o) with h^''^{s,u) determined 
by a boundary condition. For steps 2 and 4, we can then choose Ti.^ = duh^ {s,u)r"'^ sin m{9 — Oq) and T-C^ — 
/i*(s, u)r™ smm{9 — 6q). The Poisson equations in steps 5 and 6 can always be written in the form 



6±d±W, 



n-l 
p=0 



(218) 



with Wn = 0\r=R(s) ^^d 



p=0 



(219) 



with drXn = 0|j,^^^j,j, where the 0-forms Kp.n.m(s,u) and Tp,n,m{s,u) are known functions arising from previous 
iterations. Solutions to the above are 



Xr, 



^ 4(p+l)(p+l + m)' 



p=0 



2(P+1) 



m 



-R{s) 



p,n,ni{s, u) r" cosm{9 - Oq), 



2(P+1) _ ^2(p+l) 



4(p + l)(p + 1 + to) 



r smm 



(220) 
(221) 
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For example, this procedure leads to 

E{u) R'is) 



Xl_m — 



-r"rj)"cosm(6'-6lo), 

fiu) R'is) 2 _ 2. _ 

47r(l + m)i?(s)2™+i ^ ^ ^ ^0 V 

/(«) R'i-s) 



n 



47r(l + m) i?(s)2™+i 
2/(zi) i?'(s)2 



TO + 2 



47r(l + m) i?(s)2™ 
2/(w) i?'(s)2 



''•2.™ ^ 



^2 



47r(l + m) i?(s)2m 



r^r^ costo(6'-6Io), 
r'^r^ cosm{0-eo), 



fin) 



87r(l + to) V-R(s)2™+ 



R'is) 



2i?(s)2-r2 i?(,s) 



l + TO 2(2 + to) 



liO-Oo), 



R'is) 

7?(s)2™+i 
R'js)^ 2 

i?(s)2 TO 



m ^ 



TO + 2 ^/n2/2 '^ + 2- ,,2 



+ 



4(to + 2) 2to(to + 1) 



-Ris)' r 



-Ris) 



TO + 2 - 



TO 



47r(TO + 1)'" 



where f'{u) — dfiu)/du. 
Furthermore, it follows 



^'n.m = 0, 



l.m 



$2,™ =- 



2 fiu) 
m Air \Ris) 

R'is) 



r'^r'^ cosmi0-9o), 



i?'(s)2 



(222) 
(223) 
(224) 

(225) 

(226) 
(227) 
(228) 
(229) 
(230) 

(231) 
(232) 

(233) 



The multipole impedances introduced in Section [IT] are obtained by developing the above to higher order in n using 

(unsiiioi. 
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APPENDIX A: DECOMPOSITION OF 1-FORMS IN 2-DIMENSIONAL, SIMPLY CONNECTED, 

BOUNDED DOMAINS 



The method used here to analyse Maxwell's equations requires the Hodge decomposition of forms on a manifold 
with boundary (2(3| . 

Let I? be a transverse cross-section of an axially symmetric waveguide. We have the immersion map t : dV — s- V, 
so that t*a is the pull-back of a 1-form a onto the boundary dV of V. The set of smooth 1-forms on V can be 
decomposed as: 

AiiV) ^ d^TdiV) ® #^d^TiV) (Al) 

where JF(2?) is the space of smooth 0- forms on V and TdiT^) is the subspace of J^iV) whose elements satisfy the 
Dirichlet boundary condition: 

TdiV)^{f eTdiV)\L*f ^Q). (A2) 
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To demonstrate the Hodge decomposition (jAl[) we show dj^jF^(X') and ^±d±T{'D) are orthogonal with respect to 
the symmetric product: 

/ (A3) 
Jv 

where /i and ly are 1-forms on T). We will then show that if a 1-form is orthogonal to d±J-d{'D) then that 1-form must 
be in i^^d^J'iV). 

For A e J^dCD) and a G J^{V) 

(dj.A#±cUa)= / d^AA#i(#id^a) 
d±A A d^a 

i*(Adia) = (A4) 

dV 



and we conclude that d±Td{'D) and #_Ld_LJ^(X') are orthogonal with respect to (|A3I 
Let w be a 1-form orthogonal to d±!Fd{1^)- Thus, for every A e Td(V) we have 



0==(diA,cj)= / d_LAA#iW. (A5) 
However, as A satisfies the Dirichlet boundary condition, we have 



V 



0= / i*{A#^uo) 

JdV 

{d^AA#^uj + Ad^#^Lj) 

Ad^if^Lo (A6) 

and, as this is true for every A G !Fd{T>), we conclude 

di#ic^ = 0. (A7) 

The cross-section T) of the axially symmetric waveguide considered here is simply connected. By the converse of 
Poincare's Lemma (2ll |. we can therefore write 

= dj_K ^ u = -#^dj_K (A8) 

which is clearly a member of ^^d^T{T>). 
Thus, any 1-form a onV may be written as 

a = d_L^ + #±d_La (A9) 

for some A £ TdiT^) and some a E J'i'D). 



APPENDIX B: MOTIVATION BEHIND AUXILIARY POTENTIALS 



1. Maxwell Equations 

The approach used here to decompose Maxwell's equations employs the Hodge decomposition (|A9[) on transverse 
1-forms. In earlier Sections we showed that the introduction of auxiliary potentials and our approximation method 
lead to Poisson and Laplace equations for the auxiliary potentials. 

Alternatively, one may try to directly generalize the approach in Q by writing the Maxwell equations in terms of 
the electric field and postulating asymptotic series for A, a. However, as we will now show, the latter approach leads 
to third order PDEs that are harder to analyse than the hierarchy of Poisson and Laplace equations obtained earlier 
for the auxihary potentials {7^*, 7^^, 7^^ W, X}. 
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The electric field l-form E is 

E = $dC - a (Bl) 

where $ is the longitudinal component of the electric field and ~a is the transverse projection of E. The Maxwell 
system (|48ll53p can used to obtain equations for the electric field alone: 



5idi$ = 29^„$ + d^^<S>, (B3) 

^xd^a + d_L5±a = d±p + 2dLa + dha. (B4) 



2. Boundary Conditions 

The boundary of the waveguide is perfectly conducting so the electric field tangent to the boundary vanishes: 



where 



The boundary condition (IBSP with 



leads to 



S = r-R{C). (B6) 
a = ardr + aoAO (B7) 



ae = Q\r=R(Q). (B8) 
$-i?'(C)a. =0|,=H(C)- (B9) 

3. Hierarchy of approximations 

Inserting the Hodge-de Rham decomposition of a 

a = dx^ + #_Ld_La (BIO) 

into (|B2IIB4|) . where A = 0\r=R{Q, yields 

5xd±A - p - a„$ - ac$, (BU) 
(5_Ld^$ = 29^„$ + 9^^$, (B12) 
#idi<5idia + d^^^d^A = d^p + (29|„ + (d^A + #id_La) . (B13) 

As A satisfies the Dirichlet boundary condition, (IB8|) and (|B9p become 

9ra = 0|,,^fl(^), (B14) 

* + ^'(0 (Jc^ea - = 0|,=fl(^). (B15) 

Introducing an asymptotic gradual-taper expansion in the same way as in the main body of the text gives 

S^d^Ao - p, $0 = do = (B16) 

with ^0 ~ ^\r=R(s) ^"^^^ ^^o'^ " > 0, satisfies 

5^d^cE.„ - -29„|.;_i - (B17) 
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subject to the boundary condition 



r=R{s) 



Then, A„ is calculated from 

<5id^i„ = -a„i>„ - (B19) 

subject to the Dirichlet boundary condition A„ = Oj^^j^j-^-). Finally, a„ is calculated from 

#±di(5idia„ = 9„d±$„ + di$;_i + 29„dii;_i + 2a„#idid;_i + dii'„'_2 + #±dia:^_2 (B20) 

with (^rO-n = 0|r=fl(s)- Writing (|B20p in (r, 6') components yields a coupled pair of third order PDEs which is more 
difficult to tackle, in general, than the sequence of Laplace and Poisson equations for auxiliary potentials found earlier. 
However, if one can find functions and ■;/'„ such that 

di<^„ = -#±di$„, A3n = -#±dii„ (B21) 

then (jB20p can be converted to a two-dimensional Poisson equation (see equations (36, 37) in ^ for a similar approach): 

<5±dj_a„ =a„0„ + + 2a„V';_i + 2aua„_i + ?a;;„2 + K^2- (B22) 

Unfortunately, applying the transverse exterior derivative d^ to (|B2ip gives non-trivial conditions on $„ and y4„: 

= 5idi$„, = 5id_Li„. (B23) 

This implies that we can only find functions and '0„ that satisfy (|B2ip if $„ and A„ are harmonic. Thus, solutions 
to (|B21I) will not exist if the right hand sides of (jB17|) and (|B19I) are non-zero. If the right hand side of (jB19p is zero 
it follows An = since A„ satisfies Laplace's equation in TD and A„ = Oj^^^^g^. Thus, only in rare cases can we avoid 
having to solve third order PDEs. One way to avoid this problem is to employ the methods presented in the main 
text. 
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